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Boltzmann equations are often used to describe the non-equilibrium time-evolution of many-body 
systems in particle physics. Prominent examples are the computation of the baryon asymmetry of the 
universe and the evolution of the quark-gluon plasma after a relativistic heavy ion collision. However, 
Boltzmann equations are only a classical approximation of the quantum thermalization process, 
which is described by so-called Kadanoff-Baym equations. This raises the question how reliable 
Boltzmann equations are as approximations to the complete Kadanoff-Baym equations. Therefore, 
we present in this article a detailed comparison of Boltzmann and Kadanoff-Baym equations in the 
framework of a chirally invariant Yukawa-type quantum field theory including fermions and scalars. 
The obtained numerical results reveal significant differences between both types of equations. Apart 
from quantitative differences, on a qualitative level the late-time universality respected by Kadanoff- 
Baym equations is severely restricted in the case of Boltzmann equations. Furthermore, Kadanoff- 
Baym equations strongly separate the time scales between kinetic and chemical equilibration. In 
contrast to this standard Boltzmann equations cannot describe the process of quantum-chemical 
equilibration, and consequently also cannot feature the above separation of time scales. 

PACS numbers: ll.lO.Wx, 98.80.Cq, 12.38.Mh 



I. INTRODUCTION 



This article is an extension of our previous studies [H , 
where we performed a detailed comparison of Boltzmann 
and Kadanoff-Baym equations in the framework of a real 
scalar quantum field theory. The motivation for these 
studies is a better understanding of processes like lepto- 
genesis or preheating in the early universe [2, S, U [1] , 

[ or the evolution of the quark-gluon plasma after rela- 
tivistic heavy- ion collisions [1, 0, H, Q. All these phe- 
nomena require the description of many-particle systems 
out of thermal equilibrium. The standard means to deal 
with this non-equilibrium situation are Boltzmann equa- 
tions. However, it is well known that (classical) Boltz- 
mann equations suffer from several shortcomings as com- 
pared to their quantum mechanical generalizations, so- 

. called Kadanoff-Baym equations. This motivates a com- 
parison of Boltzmann and Kadanoff-Baym equations in 
order to assess the reliability of quantitative predictions 
based on standard Boltzmann techniques. In the present 
work we generalize our previous results to the case of a 
chirally invariant Yukawa-type quantum field theory cou- 
pling fermions with scalars. More precisely, we consider 
a globally SU{2) x SU{2) x U{1) symmetric quantum 
field theory, which offers two interpretations for the par- 
ticle content: On one hand, in the context of leptogenesis 
one might think of a single generation of leptons and a 
Higgs bidoublet On the other hand, this theory is 



equivalent to the linear cr-model [111, [H, [T3| , which can 
be used to describe low-energy quark-meson dynamics in 
two- flavor QCD. In any case, this work can be regarded as 
a further step to approach more realistic theories, which 
can be used to describe the phenomena motivating our 
studies. 

What are the shortcomings of Boltzmann equations? 
Originally, Boltzmann equations have been designed for 
the description of the non-equilibrium time-evolution of 
dilute gases of classical particles. As such, their range 
of validity must be scrutinized once quantum effects be- 
come relevant. This is certainly the case for elemen- 
tary particles playing the central role in phenomena like 
leptogenesis or the quark-gluon plasma. As already 
indicated above, the quantum dynamics of such sys- 
tems is described by so-called Kadanoff-Baym equations. 
Employing a sequence of approximations, Boltzmann 
equations can be derived from Kadanoff-Baym equa- 
tions [3, m, 01 [13) • However, it is important to note 
that these approximations might be neither justifiable 
nor controllable, and sometimes even inconsistent. After 
all, standard Boltzmann equations take only on-shell pro- 
cesses into account, feature spurious constants of motion, 
and introduce irreversibility by implying the ass ump tion 
of molecular chaos ( "StoBzahlansatz" ) [H, HO, [Hi HI . 
In contrast to this, Kadanoff-Baym equations are time- 
reversal invariant and take memory and off-shell effects 
into account [H, IH, HI, HI] ■ Therefore, one should per- 
form a detailed comparison of Boltzmann and Kadanoff- 
Baym equations il, iM, i25|, ,27 , 28, 29, Sfl]- 
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Due to the complexity of the problem, in a first step 
we restricted ourselves to a real scalar quantum field 
theory in 3 -)- 1 space-time dimensions fl']. Of course, in 
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FIG. 1: Closed real-time path C. This time path was invented 
by Schwinger ^331] and applied to non-equilibrium problems by 
Keldysh [SJj]. In order to avoid the doubling of the degrees of 
freedom, we use the form presented in Ref. [isj . 



this framework one can neither describe the phenomenon 
of leptogenesis nor thermalization of the quark-gluon 
plasma after a relativistic heavy ion collision. Neverthe- 
less, it certainly allowed to perform a detailed comparison 
of Boltzmann and Kadanoff-Baym equations, which re- 
vealed interesting phenomena to be investigated in more 
realistic theories. We found considerable differences in 
the results furnished by the Boltzmann and Kadanoff- 
Baym equations: On a quantitative level, we found that 
the Boltzmann equation predicts significantly larger ther- 
malization times than the corresponding Kadanoff-Baym 
equations. On a qualitative level we could verify that 
Kadanoff-Baym equations respect full late-time univer- 
sality [H, im and strongly separate the time scales be- 
tween kinetic and chemical equilibration js^ . In the case 
of a real scalar quantum field theory the Boltzmann 
equation includes only two-particle scattering processes, 
which conserve the total particle number. This spurious 
constant of motion severely constrains the evolution of 
the particle number distribution. As a result, the Boltz- 
mann equation respects only a restricted universality, 
fails to describe the process of quantum-chemical equi- 
libration, and does not separate any time scales. 

In the present work we extent our comparison of 
Boltzmann and Kadanoff-Baym equations to a chirally 
invariant Yukawa-type quantum field theory coupling 
fermions with scalars. We start from the 2PI effective 
action [sl, [s^ and derive the Kadanoff-Baym equations 
by requiring that the 2PI effective action be stationary 
with respect to variations of the complete connected two- 
point functions [s^,!!!]. First, this guarantees that the 
Kadanoff-Baym equations conserve the average energy 
density as well as global charges [s^ [13, |4l|, and sec- 
ond, the 2PI effective action has proven to be an efficient 
and reliable tool for the description of quantum fields 
out of thermal eq uilibrium in numerous previous treat- 
ments m, [H, H^, ^ . In order to derive the corre- 
sponding Boltzmann equations, subsequently one has to 
employ a first-order gradient expansion, a Wigner trans- 
formation, the Kadanoff-Baym ansatz and the quasi- 
particle approximation [11, [ll, [H, [13, [Hi • While Boltz- 
mann equations describe the time evolution of particle 
number distributions, Kadanoff-Baym equations describe 
the evolution of the complete quantum mechanical two- 
point functions of the system. However, one can de- 
fine effective particle number distributions, which can be 
obtained from the complete propagators and their time 
derivatives evaluated at equal times [H, [l^l . Finally, we 
solve the Boltzmann and the Kadanoff-Baym equations 



numerically for highly symmetric systems in 3-1-1 space- 
time dimensions and compare their predictions on the 
evolution of these systems for various initial conditions. 



II. 2PI EFFECTIVE ACTION 

We consider a globally SU{2)l x SU{2)r x U{1)b-l 
symmetric quantum field theory with one generation of 
chiral leptons and a Higgs bidoublet with Dirac- Yukawa 
interactions [l3|. The Dirac fields are denoted with 
(cc), where a is a Dirac index and / G {v^ e} denotes 
the type of the leptons. Using the Pauli matrices, the 
(complex) Higgs bidoublet can be parameterized by real 
scalar fields denoted with (x), where a G {0, . . . , 3}. 
In this notation the Lagrangian density takes the form^ 



1 



Although we refer to the scalar fields as Higgs fields, 
we would like to note that this theory is equivalent to 
the Hnear cr-model [ll|, [13, [l3l, which can be used to 
describe low-energy quark-meson dynamics in two-flavor 
QCD. This and a similar model have been considered in 
a related context in Refs. (3^. [ssj. 

Since we will compute the evolution of the two-point 
Green's functions for non-equilibrium initial conditions, 
already the classical action has to be defined on the closed 
Schwinger-Keldysh real-time contour, shown in Figure [H 
The free inverse propagators can then be read off the free 
part of the classical action: 



Iq = — J d^x d^y 
c 

1 



where the inverse free propagators are given by 

G'o.ah i^' y) = {dxi'dy^ + m%) Sc {x - y) Sab (1) 

and 

^o,Li y) ^^xSc {x ~ y) 5im ■ (2) 

We consider a system without symmetry breaking, i. e. 
($Q (x)} = 0. In this case the full connected Schwinger- 
Keldysh propagators are given by 



Gab ix,y)^{Tc{<fa (x) <^b (y)}) 



(3) 



^ The form of the kinetic term indicates that we use the Minkowski 
metric where the time-time component is negative, iti, (T2 and 
(73 are the usual Pauli matrices, while ag = il. In addition 
to the Dirac 7 matrices we will frequently use /3 = 17", and 
Pl = ^{1+ 75) and Pfl = i (1 - 75). 



2 




FIG. 2: Two-loop contribution to r2 [G,S]. Full (dashed) 
lines represent the complete connected lepton (Higgs) propa- 
gator S (G). 




FIG. 3: One-loop contribution to the Higgs self-energy. 



and 



SrIi^:y)^{Tc{^? {x)^i (y)}) 



(4) 



so that the 2PI effective action can be written in the form 

r [G, s] 



-ztrclogc [S-'] +trc [So'S] 

+ T2 [G, S] + const . (5) 

The square brackets indicate that the trace, the loga- 
rithm and the product of the propagators have to be 
taken in the functional sense, and the subscript C reminds 
us that integrals over time are running along the closed 
real-time contour. iT2 [G, S] is the sum of all two-particle 
irreducible vacuum diagrams, where internal lines repre- 
sent the complete connected propagators S and G. In 
this work we apply the loop expansion of the 2PI effective 
action up to two-loop order. The diagrams contributing 
to r2 in this approximation are shown in Figure [2l Using 
the abbreviation 



H' 



af3 
a Am 



IT] 



Q/3 



we find 



r2 [G, S] ^ -X J d^X Gaa [x, x) Gbb {x, x) 

c 

+ 2Gabix,x) Gab ix,x) 

~ 2 ^"^^ ^"'^ ^^(yUaS [x, y) HbS {y,x) ^ 



where the trace runs over Dirac and lepton type indices. 



FIG. 4: One-loop contribution to the lepton self-energy. 



They are equivalent to the corresponding self-consistent 
Schwinger-Dyson equations 



'^ab y) = ^(^oL i^^ V) - Hab {x, y) 



and 



. 5T2 [G,5] 



+ 



t^(HaS{x,y)HbS{y,x) 



and 



/ N . [G, S] 



= -KJkHt^nmSli{x,v)Gab{x,y) 
Next, we define the spectral function^ 

Ge.ab {x,y) -i([$a (a;) , $b 
and the statistical propagator"^ 

GF,ab {x, y)^l ([$a {x) , $6 (y)] + ) 



(7) 



Sim (2^' y) = iSo.ln y) ~ ^lr,i {X, y) , (8) 

where the proper self-energies are given by 



5Gba {y,x) 

= ~4XiSc {x - y) {Cdd {x, x) Sab + 2Gab {x, x) ^ 



(9) 



(10) 



III. KADANOFF-BAYM EQUATIONS 

The equations of motion for the complete propagators 
read 



dr\G,S] „ , 6T\G,S] 
— and 



SGba iy,x) 



(6) 



^ From the definition of the Higgs spectral-function we see that it 
is antisymmetric in the sense that Ggi,a{y,x) = —Ggabi^^V)- 
Furthermore, the canonical equal-time commutation relations 

give {Gg^ab{x,y))^o^yO = and {dyoG g^ah (x,y)^ = 

-SabS^ {x - y). 

^ In contrast to the spectral function, the statistical Higgs-propa- 
gator is symmetric in the sense that Gp i,a ^) = Gp ab (2^1 ?/)• 
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for the Higgs bosons, so that we can write the complete 
Higgs propagator as 

G {x, y) = Gf (x, y) - ^sign^ {x° - yO) G, {x, y) . (11) 

In the case of real scalar fields the spectral function and 
the statistical propagator are real- valued quantities [i^ . 
In a similar way, we also define the spectral function'* 



-^x^^^xJac + (x)) G g^cb {x, y) 



(12) 



and the statistical propagator^ 

Stra{^,V) = \{[n{^),^ir{y)\^) (13) 

for the leptons, so that we can decompose the complete 
lepton propagator according to 

S {x, y) = Sf (x, y) - ^sign^ - y") (x, y) . (14) 

Then, using Eqs. (fTT|) and we can decompose the 

Higgs self-energy ([9]) as well as the lepton self-energy 
according to: 

n(x,y) = (x - y) n('°™') (x) 



Hf (x, y) - -signc (x° - y°) Wg (x, y) 



and 



S [x, y) ^ T,F (x, y) - -sign^ (x° - y°) (x, y) . 

The local part of the Higgs self-energy causes a mass shift 
only, wherefore we define the effective mass by 

Af,\(x)=™y„fc + niT'^(x) (15) 

= m%Sab + 4:X(Gdd {x, x) 6ab + '2Gab (x, x) ^ . 

After convoluting Eqs. ([7|) and ([8]) from the right with 
the corresponding complete propagators, we observe that 
both equations split into two complementary evolution 
equations for the statistical propagators and the spectral 
functions, respectively 38]: 

{-dxi-dx^6ac + (x)) GF,cb {x, y) 
d'^z TlF.ac {x, z) Gg^cb {z, y) 



- j d'^zUg^ac{x,z)GF,cb{z,y) , (16) 





The adjoint of the lepton spectral-function is given 
•5'e,im(^'f) = -PSe,ml{y,x) [3. Furthermore, the 
canonical equal-time anti-commutation relations give 

^ The adjoint of the statistical lepton-propagator is given by 

sii,„{x,y) = f^SF,rni {y,x)i3. 



d'^zIlg,acix,z)Gg,cb{z,y) , (17) 



vSFJm {X, y) ^ J '^FJk (X, z) Sg^km {z , y) 





d'^z T.g^ik (x, z) SF,km (z, y) 



(18) 



and 



cSg^im {x, y) 



d'^zT.g^ik{x,z)Sg^k7n{z,y) ■ (19) 



Nowadays, it is practically impossible to solve the Ka- 
danoff-Baym equations numerically in this general form. 
However, for initial conditions which are invariant un- 
der spatial translations, spatial rotations, parity, charge 
conjugation, and chiral transformations, the propagators 
take the form [33 



G, 
and 



ab{x,y)^ /(^~j^ exp(^ifc(a;-y)^G(x°,y°,fc)^a6 

(20) 



Sim {x, y) = / -^-^ exp (ik {x ~ y) 
(27r) ^ 



(21) 



X (^S'i.{x°,y'',k)jo-i^Sv{x",y'',k)^j^ 5im , 

where k = |fc|. The index V indicates that Sy would 
transform as a vector under a Lorentz transformation. 
Due to CP invariance, the statistical and spectral vector 
components of the lepton propagator satisfy 

•^y.F (■^^' 2/*^' ^) ^ ^ Sy p {jj'^ T x^ , k^ = Sy pl^x^jy^jk"^ 

Sv,F {x'>,y",k) = Sv.F {y'\x",k) = S'^^p {x°,y'>,k) 

50,,(xO,y",A:) = 50,JyO,x",fc) = 5£ (x", y". A:) 

Sv,g{x°,y°,k) = ~Sv.g{y°,x',k) = S*y ^ {x° , y° , k) 

(22) 

Thus, the explicit factor of i makes Sy Sv,f, Sy ^ and 
Sv,g real- valued quantities: 

Im {S^p) = Im {Sv,f) = Im {S°y g) = Im {Sv,g) - . 

(23) 

Furthermore, the canonical equal-time anti-commutation 
relations for fermion fields imply 

S^„ (t, t,k) = l and Sv.g (t, t,k) = 0, (24) 
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the second equality being consistent with the anti- 
symmetry of the space-Uke vector-component of the lep- 
ton spectral function, cf. Eq. ((22|) . Of course, the rela- 
tions (HOI), (EH, (1221) and ^ also hold for the corre- 



sponding self energies, so that the Kadanoff-Baym equa- 
tions can be simplified drastically. The simplified Kada- 
noff-Baym equations for the Higgs propagator read [sHi 



y" 

= J dz° Uf {x°, k) Gg (z°, y°, k)- j dz° Ilg {x°, z°, k) Gp k) 





(25) 



and 



[dlo +e + (x°)] G, yO, k)=- dz^ n, {x\z\ k) G, y\ k) 



(26) 



In the same way, the 128 complex-valued Kadanoff-Baym equations (US]) and for the lepton propagator can be 
reduced to the following 4 real- valued equations fssj : 



dz° 



dz'' 



i:°y^F {x°, z\ k) S^y^^ y\ k) + ^v,F z°, k) Sv,g y\ k) 
(x", zO, k) S^,F (z", A:) Ey,, (xO, z°, k) Sv,f {z°, y°, k) 



(27) 



Ev,^ z", k) (zO, y", fc) - S°._^ z°, A;) Sv,e {z°, y° , k) 



V 







dz^ 



Ev,, (x", z°, A:) 5^^,^ (z", k) - E?.,, (x", z", A:) Sv,f (z°, y", fc) 



(28) 



and 



d^oSl,^{x\y\k) + kSv..s{x\y\k) 



?o /^O „,0 

■Q 
x° 



EO,,^ z°, k) (zO, y", k) + Ey,, z", k) Sv,e {z\ y\ k) 



>;0 



5.o5y,,(xO,y°,A;) -fc^o.,, (x°,y°,A:) 



dz" 



Ey, {x\ z", fc) 5°,, (z", y", fc) - E?.,, (x°, z°, fc) (z", y", k) 



(29) 



(30) 
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The expressions for the Higgs and lepton self-energies are 
given in App.El As explained in more detail in Ref. [s^, 
one can define an effective kinetic energy distribution 
uj{t,k), as well as effective scalar and fermion particle 
number distributions Ug {t, k) and nf (t, fc), which can be 
obtained from the statistical propagators according to 



Gf ix°,y°,k) 



and 



Us {t,k) ^uj{t,k) Gf it,t, k) 



Uf {t, ^) = 2 ^ "^^'^ 



(31) 



(32) 



(33) 



The definition of such particle numbers is necessary in 
order to make contact to Boltzmann equations, e. g. when 
comparing numerical solutions of Boltzmann and Kada- 
noff-Baym equations, which we wiU do in Sect. El We 
emphasize, however, that the Kadanoff-Baym equations 
are self-consistent evolution equations for the complete 
propagators of our system, and that one has to follow 
the evolution of the two-point functions throughout the 
complete a;°-?/°-plane (of course, constrained to the part 
with a;° > and > 0). One can then follow the 
evolution of the effective particle number densities along 
the bisecting line of this plane. 



IV. BOLTZMANN EQUATIONS 

In this section we briefly sketch the standard way of de- 
riving Boltzmann equations from Kadanoff-Baym equa- 
tions jli,[il[il,|4l|4i,|43. One has to employ a Wigner 
transformation, a first-order gradient expansion, the Ka- 
danoff-Baym ansatz and the quasi-particle approxima- 
tion. 

First, we subtract Eq. ^ (Eq. ^) with x" and 
y^ interchanged from Eq. (^5)1 (Eq. ^7}). Then we re- 
parameterize the propagators and the self energies by 
center and relative times, e. g. 



G{u°,v°,k) =G 



Next, we define the center time t = ^ and the relative 
time = a;° — and observe on the left hand side of 
the difference equations that 

9„o9t.o — d,,od,,o = 2dfd^o 



and 



dyO 



dt 



are automatically of first order in dt- Furthermore, we 
Taylor expand the effective masses on the left hand side 



of the difference equation for the scalars as well as the 
propagators and self energies on the right hand sides of 
both difference equations to first order in dt around t. 
After that, we Fourier transform the difference equations 
with respect to the relative time s°. The Wigner trans- 
formed scalar statistical propagator and scalar spectral 
function are given by 



GF{t,uj,k) = J ds" e^p{iujs")GF {t,s",k) , 
Gg{t,Lo,k) ^ ~i [ ds° exp{iujs°)Ge{t,s°,k) 



As Gg (t, s°,/c) is a real- valued odd function of the rel- 
ative time s*^, we introduced an explicit factor of —i in 
order to make its Wigner transform again a real-valued 
function. For similar reasons we also introduce a fac- 
tor of —i for the Wigner transforms of Sv,g {t, , k) , 
S% {t, s°,k), S°. ji (i, s°,k), and 8°,,^ {t, s", fc)^, as well 
as the corresponding self energies. In order to be able 
to really perform the Fourier transformation, we have to 
send the initial time to — oo. At least for large x^ and 
this can be justified by taking into account that cor- 
relations between earlier and later times are suppressed 
exponentially ^, [42]. For early times, however, this is 
certainly not the case. The result of all these transfor- 
mations are quantum-kinetic equations for the statistical 
propagators Gf and S^p [16, 17, 18, 48, 49, 50, 51]^: 



{O; Gf} 



PB 



TloGi 



liFGg + iliF] Re(Gfl,)} 



PB 

(34) 



and 



{W; Syp^pg — Yiy^^Syp — TfypSy^ + Tiy^gSy^F 

- ^v,fSv^s - Re {S%)} (35) 

+ {Re (Sy^fl) ; Sy^F}pB + {^V,F', Re(5y^fl)}p^ , 

where the Poisson brackets are defined by 

if'^ 9}pB ^ - 9tf d^g + d^f dtg . 
The auxiliary functions 

n {t, Lj, k) = -uj^ + fc2 + (t) + Re {Ur {t, uj, k)) 
and 

W (t, w, fc) = w + Re {T."y n (i, w, k)) 



® The retarded and advanced propagators, e. g. Gn (x" , j/" , fc) = 
e{x'' -y°)Ge{x°,y'^,k) and GA{x°,y°,k) = 
—9[y^—x^)Gg(^x'^,y'^,k), and self energies have to be 
introduced in order to remove the upper boundaries of the 
memory integrals in the Kadanoff-Baym equations. 

^ The complete and closed set of these quantum-kinetic equations 
comprehends 9 equations and self energies, which, for complete- 
ness, are shown in App. iBl 
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have been introduced to simplify the notation. Employ- 
ing the first-order Taylor expansion is clearly not jus- 
tifiable for early times when the equal-time propagator 
is rapidly oscillating [l|, [2^. Consequently, one might 
expect that the above quantum-kinetic equations and 
also the Boltzmann equations, which we derive subse- 
quently, fail to describe the early-time evolution and that 
errors accumulated for early times cannot be remedied at 
late times. In fact, the first-order gradient expansion is 
motivated by equilibrium considerations: In equilibrium 
the propagator depends on the relative coordinates only. 
There is no dependence on the center coordinates, and 
one may hope that there are situations where the propa- 
gator depends only moderately on the center coordinates. 
This is clearly the case for late times when our system is 
sufficiently close to equilibrium. However, already after 
moderate times the rapid oscillations mentioned above, 
have died out and are followed by a monotonous drifting 
regime P, H^l . In this drifting regime the second deriva- 
tive with respect to t should be negligible as compared 
to the first-order derivative, and a consistent Taylor ex- 
pansion can be justified even though the system may still 
be far from equilibrium. Here, it is crucial that the Tay- 
lor expansion is performed consistently for two reasons: 
First, this guarantees that the quantum-kinetic equations 
satisfy exactly the same conservation laws as the full 
Kadanoff-Baym equations do [l7[. Second, it has been 



shown that neglecting the Poisson brackets severely re- 
stricts the range of validity of the quantum-kinetic trans- 
port equations 51, 52]. 

In order to derive Boltzmann equations from the 
quantum-kinetic equations for the statistical propagators 
([M)) and ([55)) . first we have to discard the Poisson brack- 
ets on the right-hand sides, thereby sacrificing the consis- 
tency of the gradient expansion. On the left-hand sides 
we remove the time dependence of the auxiliary quanti- 
ties fl and W. We take 



and 



W {t,uj,k) =LU 



n{t,uj,k) = -uj"^ + k'^ + ' 



where m is the thermal mass of the scalars. After that, 
we employ the Kadanoff-Baym ansatz 

Gf (t, uj, k) = Gg (w, k) Q + Us (t, uj, fc) j (36) 



and 



Sf {t, uj, k) = Sg {uj, k) Q - 71/ (t, u, fc)^ 



(37) 



which also can be motivated by equilibrium considera- 
tions. In fact, this is a generalization of the fluctuation- 
dissipation theorem, which states that, for a system in 
thermal equilibrium, the statistical propagator is propor- 
tional to the spectral function. The fluctuation dissipa- 
tion theorem can be recovered from Eqs. and ([77)1 



by discarding the dependence on the center time t and 
fixing Ug and riy to be the Bose- Einstein and Fermi-Dirac 
distribution function, respectively. The last approxima- 
tion, which is necessary to arrive at the Boltzmann equa- 
tions, is the so-called quasi-particlc (or on-shell) approx- 
imation. For the scalars this means that the spectral 
function takes the form 

Gg{uj,k)^-^^{5{u:-E{k))~5{Lo + E{k))) , 
where the quasi-particle energy is given by 



E{k) = Vm^ + P _ 
For the lepton spectral function we assume 

S^^^ (w, k) = 7r{6iLj-k) + S{uj + k)) , 
Sv,gi(^,k) ~ TT (S {ll> — k) — S (uj + k)) . 

Once more, we would like to stress that the exact time 
evolution of the spectral functions is determined by the 
Kadanoff-Baym equations. It has been shown that the 
spectral function can be parameterized by a Breit-Wigner 
function with a non- vanishing width [261 . [29| . To reduce 
the width of this Breit-Wigner curve to zero is certainly 
not a controllable approximation and leads to very large 
qualitative discrepancies between the results produced by 
Kadanoff-Baym and Boltzmann equations. In fact this 
approximation can only be justified if our system con- 
sists of stable, or at least very long-lived, quasi-particles, 
whose mass is much larger than their decay width. We 
also would like to note that a completely self-consistent 
determination of the thermal mass in the framework of 
the Boltzmann equation requires the solution of an inte- 
gral equation for E (k), which would drastically increase 
the complexity of our numerics. As none of our physical 
results depend on the exact value of the thermal mass, for 
convenience, we use the equilibrium value of the thermal 
scalar mass as determined by the Kadanoff-Baym equa- 
tions. Eventually, we define the quasi-particle number 
densities by 



and 



Us (t, k) = Us {t, E (k) , fc) 



nf {t,k) = Uf (t, k, k) 



After equating the positive energy components in 
Eqs. (|34p and ((35)l we arrive at the following Boltzmann 
equations: 

dtn, (t, k) = 2nrf [ / d^q S^k - p - q) 

J (2it) J 



X S{E{k) -p- q)' 



pq 

pq 



m 

(^Us (<, k) + ij Hf (t, p) nf {t, q) 
■s [t, k) {rif {t,p) - {rif {t, q) 



(38) 



7 



and 



dtUf [t, k) = I / d\ 5^{k+p-q) 

J (zTTj J 



x5{k+p-E{q)) 



1 A fcp 

S(q) kp 



(jif (t, k) - 1^ (uf {t,p) - l^ns (t, 9) 
n/ (i, fc) nj {t,p) (us {t, q) + 1^ 



(39) 



Exploiting isotropy the above 6-dimensional Boltzmann 
collision integrals can be reduced to 1-dimensional inte- 
grals: 



dtHs (t, k) = 



dp (go) Js{k,p, go) 



Tr^E{k) 


Us [t, k) + ijuf {t,p) Uf {t, go) 
{t,k) (uf (t,p)~l) {n.f (t, (7o)"l 



(40) 



2rf f g 
dtUf {t, fc) = — J dgd (pq) -^J^Jf ik,Po, g) 



(uf {t, k) - 1^ {nj {t,po) - ijn^ (t, g) 



Uf (t, k) Uf {t,po) (us {t, g) + 



(41) 



The details of these calculations and the definitions of all 
of the auxiliary quantities are given in App.O This sim- 
plification of the collision integrals is crucial in order to 
implement efficient computer programs for the numerical 
solution of the Boltzmann equations. 

In this section we showed that the derivation of Boltz- 
mann equations from Kadanoff-Baym equations requires 
a number of non-trivial approximations and assumptions. 
One has to employ a first-order gradient expansion, a 
Wigner transformation and a quasi-particle approxima- 
tion. In this sense, one can consider the Kadanoff-Baym 
equations as quantum Boltzmann equations re-summing 
the gradient expansion up to infinite order and including 
memory and off-shell effects. 



V. COMPARING BOLTZMANN 
VS. KADANOFF-BAYM 

A. Initial Conditions and Numerical Settings 

In order to solve the Kadanoff-Baym equations numer- 
ically, we follow exactly the lines of Refs. [l|, IH, on 
a lattice with 2000^ x 32^ lattice sites. The values for 



the coupling constants are 77 = 1 and A = 0.25. The 
initial conditions for the statistical propagators are de- 
termined by scalar and fermionic particle number distri- 
butions Hs (t ~ 0,p) and Uf {t = 0,p) according to 



Gf {x°,y°,p)^^^ 



Us (i = 0,p) + i 



^^0 {p) 



[d.oGF{x°,y",p)].,^^ = 



= ^0 {p} (ns {t ^0,p)+ ^ 



2 ' 



^V.F i,^^ ^iP ^P) x°=y°=0 ~ 

Sv,F {x°,y°,p)^o^y0^o = ^-nf{t = 0,p) , 

where loq (p) is the initial scalar kinetic energy distribu- 
tion. On the other hand, the initial conditions for the 
spectral functions are determined by equal-time (anti-) 
commutation relations. We solve the Boltzmann and Ka- 
danoff-Baym equations for three different sets of initial 
particle number distributions, which are shown in Fig- 
ure [5l All initial conditions correspond to the same (con- 
served) average energy density. Above that, for the initial 
conditions ICl and IC2 also the sum of the initial scalar 
and fermionic average particle number densities agree. 
The numerical solution of the Boltzmann equations pro- 
ceeds along the lines of Ref. 1]. The scale in our plots 
is set by the scalar thermal mass m = (p = 0), where 
ujeq (p) is the effective kinetic energy distribution pTjl for 
sufficiently late time t. 



B. Universality 

Figs. [5] and [7] exhibit that the Kadanoff-Baym equa- 
tions respect full universality: Figure [S] shows the time 
evolution of the particle number distributions for a fixed 
momentum mode. The particle number distributions 
start from different initial values and go through very 
different early-time evolutions. Nevertheless, in the case 
of the Kadanoff-Baym equations they all approach the 
same universal late-time value. Figure [7| shows the parti- 
cle number distributions for times when equilibrium has 
effectively been reached. In the case of Kadanoff-Baym 
equations, we observe that the equilibrium number dis- 
tributions agree exactly independent of the initial condi- 
tions, which proves that we could have shown the plots 
of Figure [5] for any momentum mode. In particular, the 
straight lines in Figure [7] prove that the equilibrium num- 
ber distributions take the form of Bose-Einstein or Fermi- 
Dirac distribution functions with a universal temperature 
T = 2.3m and universally vanishing chemical potentials. 

In contrast to this, Boltzmann equations maintain only 
a restricted universality. Figs. [S] and [7] reveal that only 
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— ICl 

— IC2 
ICS 



0123456 0123456 

p I m p I m 

FIG. 5: Initial particle number distributions for scalars (left) and fermions (right). These distributions define the initial 
conditions (IC) for which we numerically solved the Boltzmann and Kadanoff-Baym equations. All initial conditions correspond 
to the same (conserved) average energy density. Above that, for the initial conditions ICl and IC2 also the sum of the initial 
average number densities agree. 



the initial conditions ICl and IC2 lead to the same late- 
time behavior, which deviates significantly from the one 
approached by the third initial condition ICS. Again, 
the straightness of the lines in Figure [7] proves that the 
equilibrium number distributions take the Form of Bose- 
Einstein or Fermi-Dirac distribution functions. However, 
the different slopes of the lines indicate that the temper- 
ature is not the same for all initial conditions, and the 
non-vanishing y-axis intercepts indicate that Boltzmann 
equations may predict different non-vanishing chemical 
potentials. Fitted values for these quantities are given in 
Table H 

The reason for the observed restriction of universality 
can be extracted from Figure [8] where we plotted the 
time evolution of the average particle number densities 
per degree of freedom 



and 



Nfit) 



(2^ 



(42) 



(43) 



and their sum. Provided <I>-derivable approximations are 
employed, Kadanoff-Baym equations conserve the aver- 
age energy density as well as global charges [H, HO, EH ■ 
However, as we consider systems with vanishing net 
charge density neither of the above average particle num- 
ber densities has to be conserved, nor their sum. Indeed, 
Kadanoff-Baym equations include off-shell particle cre- 
ation and annihilation [l], [2^ , so that all of the quanti- 
ties plotted in Figure [8] may change as time goes on and 
approach a universal equilibrium value. 

In contrast to this, due to the quasi-particle (or on- 
shell) approximation the Boltzmann equations (|40p and 
(I4ip only include decay and recombination processes of 





T,/m 


Tf/m 


fj.a/m 






ICl 


1.044 


1.044 


0.910 


0.455 


1.999 


IC2 


1.044 


1.044 


0.910 


0.455 


2.001 


ICS 


1.125 


1.124 


0.586 


0.293 


1.997 



TABLE I: Temperatures and chemical potentials as predicted 
by the Boltzmann equations. The values in this table have 
been obtained by fitting the equilibrium particle number dis- 
tributions shown in Figure [7] against Bose-Einstein and Fermi- 
Dirac distribution functions, respectively. 



the form 



1 scalar 



2 fermions 



(44) 



More precisely, one of four scalars may decay into or be 
recombined from one of two fermion pairs. As a con- 
sequence the sum of the average particle number densi- 
ties p2|) and ([43]) is strictly conserved, as can be seen 
in Figure [51 Of course, this artificial constant of mo- 
tion severely restricts the evolution of the particle num- 
ber distributions. As a result, the Boltzmann equations 
maintain only a restricted universality and, as will be dis- 
cussed in the next subsection, fail to describe the process 
of quantum-chemical equilibration. 



C. Chemical Equilibration 

In a system allowing for creation and annihilation of 
particles, the chemical potential of particles, whose to- 
tal number is not restricted by any conserved quantity, 
must vanish in thermodynamic equilibrium. Accordingly, 
as we consider systems with vanishing net charge density 
the chemical potentials for scalars and fermions should 
vanish once equilibrium has been reached. Indeed, Kada- 
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FIG. 6: Time evolution of the particle number distributions for a fixed momentum mode p. We see that Kadanoff-Baym 
equations respect full universality, whereas Boltzmann equations maintain only a restricted universality. 



noff-Baym equations lead to universally vanishing chem- 
ical potentials. In contrast to this, as one can see in 
Table IH the Boltzmann equations in general will give 
non-vanishing chemical potentials. For a system which 
includes only interactions of the form in equilib- 

rium the chemical potentials are expected to satisfy the 
relation 

As one can see in the right-most column of Table HI this 
relation is indeed fulfilled up to numerical errors < 0.3%. 
Thus, the Boltzmann equations and (|4ip lead to 
a classical chemical equilibrium. As mentioned above, 
however, quantum-chemical equilibrium requires that the 
chemical potentials vanish for systems with vanishing net 
charge density. In this sense, the non-vanishing chem- 
ical potentials in Table U indicate that the description 
of quantum-chemical equilibration is out of reach of the 
Boltzmann equations ([^0]) and (|¥T|) . 

D. Separation of Time Scales 

As has been discussed in Ref. [1[ in the framework 
of a purely scalar theory and in Ref. (3^ in the frame- 



work of the linear sigma model underlying our studies 
in this work, Kadanoff-Baym equations strongly sepa- 
rate the time scales between the kinetic and the complete 
thermodynamic (including chemical) equilibration. This 
phenomenon has been called prethermalization [32] , and 
implies that certain quantities approach their equilibrium 
values on time scales which are dramatically shorter than 
the thermodynamic equilibration time. 

As we have seen in this work and in Ref. [ij, stan- 
dard Boltzmann equations cannot describe the phe- 
nomenon of quantum-chemical equilibration, and thus 
they also cannot describe the approach to the quantum- 
thermodynamic equilibrium. Consequently, standard 
Boltzmann equations cannot separate the time scales be- 
tween the kinetic and the full thermodynamic equilibra- 
tion and hence a description of prethermalization is out 
of reach of standard Boltzmann equations. 

VI. CONCLUSIONS 

In this article we addressed the question how reliable 
Boltzmann equations are as approximations to Kadanoff- 
Baym equations in the framework of a chirally invariant 
Yukawa-type quantum field theory coupling scalars with 
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FIG. 7: Equilibrium particle number distributions. The straight hnes in these plots prove that in equilibrium the particle 
number distributions indeed take the form of Bose-Einstein or Fermi-Dirac distribution functions, where the temperature is 
given by the inverse slope and the chemical potential is proportional to the y-axis intercept. The Kadanoff-Baym equations lead 
to a universal temperature T = 2.3m and universally vanishing chemical potentials. The temperatures and chemical potentials 
predicted by the Boltzmann equations are given in Table U 



fermions. Starting from the 2PI effective action, we re- 
viewed the derivation of the Kadanoff-Baym equations 
and the approximations which are necessary to eventu- 
ally arrive at standard Boltzmann equations. We solved 
the Boltzmann and Kadanoff-Baym equations numeri- 
cally for highly symmetric systems in 3-1-1 space-time di- 
mensions without any further approximations and com- 
pared their solutions for various non-equilibrium initial 
conditions. 

We demonstrated that the Kadanoff-Baym equations 
respect universality: For systems with equal average en- 
ergy density the late-time behavior coincides indepen- 
dent of the details of the initial conditions. In par- 
ticular, independent of the initial conditions the parti- 
cle number distributions, temperatures, chemical poten- 
tials and thermal masses predicted for times, when equi- 
librium has effectively been reached, coincide. Above 
that, Kadanoff-Baym equations incorporate the process 
of quantum-chemical equilibration: For systems with 
vanishing net charge density the chemical potentials van- 
ish once equilibrium has effectively been reached. Last 
but not least, Kadanoff-Baym equations feature the phe- 
nomenon of prethermalization and separate the time 
scales between kinetic and full thermodynamic (including 



quantum-chemical) equilibration. 

The quasi-particle approximation introduces spurious 
constants of motion for standard Boltzmann equations 
(cf . Figure [5]) , which severely restricts the evolution of 
the particle number distributions. As a result, Boltz- 
mann equations cannot lead to a universal quantum- 
thermal equilibrium and maintain only a restricted uni- 
versality: Only initial conditions for which the average 
energy density, all global charges and all spurious con- 
stants of motion agree from the very beginning, lead 
to the same equilibrium results. As shown in Table [H 
Boltzmann equations cannot describe the phenomenon 
of quantum-chemical equilibration and, in general, will 
lead to non-vanishing chemical potentials even for sys- 
tems with vanishing net charge density. Due to the lack 
of quantum-chemical equilibration, the separation of time 
scales observed for the Kadanoff-Baym equations is ab- 
sent in the case of Boltzmann equations, which renders 
the description of prethermalization impossible. 

Some of the approximations, which are required to 
derive Boltzmann equations from Kadanoff-Baym equa- 
tions, are clearly motivated by equilibrium considera- 
tions. Taking the observed restriction of universality 
into account, we conclude that in the context of relativis- 
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FIG. 8: Time evolution of the average particle number densities and their sum. KadanofF-Baym equations allow for a change of 

the latter quantity, whereas it is artificially conserved by Boltzmann equations. The quantitative disagreement of the average 
particle number densities can be attributed to different initial conditions and discretization schemes for the Boltzmann and the 
Kadanoff-Baym equations, and is of no relevance for the purpose of the present work. 



tic quantum fields one can safely apply standard Boltz- 
mann equations only to systems which are sufficiently 
close to equilibrium, so that the spurious constants of 
motion emerging in Boltzmann equations already take 
their equilibrium values. However, for systems far from 
equilibrium standard Boltzmann equations work reliably 
neither for early times (no prethermalization) nor for late 
times (only restricted late-time universality, no quantum- 
chemical equilibration). Accordingly, for systems in the 



intermediate regime the results given by standard Boltz- 
mann equations should be treated with care. For realis- 
tic scenarios, like leptogenesis or the quark- gluon plasma, 
non-negligible corrections to Boltzmann equations are ex- 
pected, which should be evaluated. 

Solving Kadanoff-Baym equations numerically is sig- 
nificantly more difficult than solving the corresponding 

standard Boltzmann equations. However, the consider- 
able discrepancies found for numerical solutions of Ka- 
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danoff-Baym and Boltzmann equations revealed equally 
significant limitations for standard Boltzmann equations. 
Accordingly, the importance of numerical solutions of 
Kadanoff-Baym equations cannot be over-estimated and 
it is certainly worth to face the arising difficulties. 

In the present work we considered standard Boltzmann 
equations at lowest order in the particle number densities, 
and we employed the standard Kadanoff-Baym ansatz 
for their derivation. Further studies are needed in or- 
der to estimate whether and in how far the situation for 
Boltzmann equations can be improved by including non- 
minimal collision terms or by employing a generalized 
Kadanoff-Baym ansatz. 

In the future it will be important to perform a similar 
comparison of Boltzmann and Kadanoff-Baym equations 
also in the framework of gauge theories. Above that a 
complete quantum mechanical description of leptogene- 
sis would require a treatment of Kadanoff-Baym equa- 
tions on an expanding space-time, which induces fur- 
ther non-equilibrium effects. Independent of the com- 
parison of Boltzmann and Kadanoff-Baym equations 
we are looking forward to learn to which extend an 
entirely non-perturbative renormalization procedure af- 
fects the results quantitatively. Above all, such a non- 
perturbative renormalization procedure should have a 
stabilizing virtue for the computational algorithms. 
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APPENDIX A: SELF ENERGIES 

In this appendix we give the expressions for the self en- 
ergies, which have to be inserted in the simplified Kada- 
noff-Baym equations to According to Eq. 
the effective mass in Eqs. (pSj) and (^5]) is given by 

- ml + 24A j -^Gp {x\x\p) . 

Using the notation 

the statistical and spectral Higgs self-energies can be 
written in the form 



and 



k) = -8,7^ j 1p f d\5{k-p-q) 

J (ZTT) J 



Sv,F y°, q) Sv,F - ^Sv,, y', q) Sv,^ {x", y°,p) 



n, {x°, y\ k) = -mf J^Jd'qS{k-p-q) 

- 5°,, (x", y", q) S% {x', y\p) + Sv„ {x^ , y\ q) Sv,F (x", 



The simplified lepton self-energies are given by 



y\ k) ^ -4^^ [ 1p [ d^q5{k~p~q) 

J (ZTT) J 



Gf {x',y\q) S^y^p {x'^y\p) - \g, {x\y^q) 5°,, {x\y\p) 
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El^^ k) = -V j ^ j d\5{k-p-q) 



and 



APPENDIX B: QUANTUM-KINETIC 
EQUATIONS 

Here we give the complete set of quantum kinetic equa- 
tions, which are obtained from the simplified Kadanoff- 
Baym equations to ([50)1 once one performs a Wigner 
transformation and a first-order gradient expansion. The 
quantum-kinetic equations for the scalars read: 

- {ri; Gf}pb = ^qGf - ^fGq + {Uf; Re(Gi^)}p^ , 



-{n; G,}pg = {n,; Re{GR)}p^ 



Gf 



(Bl) 



— {Re (Sy^ij) ; 5'y,e}p^ - {^v.e'i R'^iSv,R)}pB ' 



and 



{W; Svjpg = {j:"y^^; Re{Sv,R)} 



PB 



{Re {Ev..r) ; ^y,,}ps + {^v^.^; Re (5°,^) } 



PB 



The quantum-kinetic equations for the retarded lepton 
propagators read 



cO _ 



{W - \YPy^^' - (fc - Re (Eka.) - f Sy,,) 



It can be shown that 116 



= 21m(Gp) = 



^2 



in2 

i g 



indeed satisfies the kinetic equation for the scalar spectral 
function (jBip . The quantum-kinetic equations for the 
fermionic statistical propagators and spectral functions 
read: 



and 



fc - Re - f Sv,g 

{W - iEO.^' - (fc - Re - ^T^v.,) 



^V,F'-'V,Q 



V,Q^V,F 



- T,v.FSv,g ~ {Sy Re (5*^ p) } 

+ {Re(i;yfl); S'ypjp^ + {Syp; Re (5'y,fl)}p^ , 



APPENDIX C: SIMPLIFYING THE 
BOLTZMANN COLLISION INTEGRALS 



{W; Sv,f}pb = ^v.Q^v.F + ^v,F^v.Q 



S?/ p'S'y.o + ^V.Q^v.F 



^V.pSv.g - 



PB 



{S0.,p; Re{SvM))pB + {^y^F-^ Re(5°^p)} 



PB 



This appendix reveals the details of the calculation 
leading from the Boltzmann equations ([38]) and ([39|) 
to their simphfied versions and (PT|) |ll. [53 |. For 
zero momentum the evaluation of the collision integral 
in Eq. ([55]) is literally trivial: 
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dtn,{t,k = 0) = - 



mrj 



(us (t, k = 0) + l^nf [t, p) Uf (t, q) 

(t, fc = 0) (uf {t,p) - l) (uf {t, q) - l) 



For fc > a little more work has to be done. We rewrite 
Eq. ([55]) using the Fourier representation of the momen- 
tum conservation 5 function 

6^{k-p-q)^ / -—3 exp + ip^ + iq^) 

J (2n) 

and spherical coordinates. The scalar product of two 
vectors is then given by 

pq = pq I sin t9p sin i)q cos (fp — ipq) + cos cos dq j . 



We perform the integrals over the solid angles in the order 
Clq, flp, fl^. Using the notation 



■ / ^ sin {x) 

] (X) — cos (x) 

X 



we find 



J dVlq exp (iqi) - 1^ = ^[i cos {-dp) j {q£_) - sin {q£_) 
j dQ.p exp (ip^) (^i cos {-dp) j (q^) - sin (gO ) = ^P^'^ ^ ('?^) + ^P^^ ^^"^ ('?^) 



47r 

dfi^ exp (— jfc^) = — sin(fc^) . 



After defining the auxiliary function 



Js (fc, p, q) - ~pq I ^^^^ (j (K) j {qO + Sin {pO sin {qO ) 







lifc(^^ - (P + 9)^ ) (sign (fc - p - g) - sign (k + p - q) 
- sign (fc - p + g) + sign (fc + p + g) ^ 



and integrating over g, we eventually arrive at Eq. ()40[) . 
where 

go = £^ (fc) - P • 

Next, we work out the collision integral for the fermions. 
First of all, we integrate Eq. (j39|) over il^. On the left 



hand side this gives a factor of 47r. On the right hand 
side we evaluate the integrals over the solid angles in the 
order 17^, fip, Ofc, O^: 
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f An 

J dilq exp(zql) = — sin (q^) , 

J d% exp {-ipC) + 1^ = ^ ( sin (P^) ~ * (^fe) 'yP^ 

j d^k exp {-iki) (sin (p^) - i cos (i^fc) j (p^) = ^ ( (^^) '^^'^ (f^^) ~ (^^) (p'^) ) ' 



dQ.e = 47r 



Defining the auxiliary function 



Jf{k,p,q) ^ p j di ( sin (fcO sin (K) - {kS) J (K) ) 




TT 



16A:2 ( '•'^ ^ ^-"^ ^ '^^) ('^^^'^ (k-p-q) - sign (fc + p - g) 
- sign (fc — p + g) + sign (fc + p + g) ^ 

I 



and integrating over p yields Eq. (I41I1 . where 
Po^ E(q)~k . 
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